Gauge (non-) invariant Green functions of Dirac fermions coupled to gauge fields 



OO 



O 



O 

o 



D. V. Khvcshchenko 

Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599 

We develop a unified approach to both infrared and ultraviolet asymptotics of the fermion 
Green functions in the condensed matter systems that allow for an effective description in the 
framework of the Quantum Electrodynamics. By applying a path integral representation to the 
previously suggested form of the physical electron propagator we demonstrate that in the massless 
case this gauge invariant function features a "stronger-than-a-pole" branch-cut singularity instead 
of the conjectured Luttinger-like behavior. The obtained results alert one to the possibility that 
construction of physically relevant amplitudes in the effective gauge theories might prove more 
complex than previously thought. 
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In a generic many-body fermion system, a repulsive electron-electron interaction is normally expected to result 
, in a suppression of any amplitude which describes propagation of fermionic quasiparticles. For instance, in the 
phenomenological Fermi liquid theory, the residue of the electron Green function G(e,p) = Z(e)/(e — E(p) + fi) gets 
reduced compared to the non-interacting value (Z(e) = 1), thus exhibiting a partial (0 < Z(0) < 1) suppression of 
the simple pole which corresponds to the bare fermionic quasiparticles. 

The question as to whether or not the repulsive fermion interactions can result in an even more severe, complete, 
destruction of the pole (Z(0) = 0) remains a subject of an ongoing debate. Such a behavior is well known to occur in 
the one-dimensional (ID) Luttinger and related models with short-ranged interactions, in which case the residue of 
the fermion Green function exhibits a characteristic algebraic decay Z(p) ~ p 11 as a function of the Lorentz- invariant 
£SJ , momentum p = y— p 2 — (p 2 — uj 2 ) 1 / 2 and is controlled by an anomalous dimension r\ > 0. 

In the ID coordinate space, this behavior corresponds to the suppression of the electron propagator G(t, x) ~ 
exp(±i&Fx)/|x ± t\ 1+ri which at long times and distances decays faster than the non-interacting one (77 = 0). In 
the absence of spin, the above Green function is Lorentz invariant, apart from the oscillating factors exp(±ikpx) that 
stem from a finite (2k p) separation between the two ID Fermi points, in accordance with the fact that the low-energy 
excitations iPr,l confined to the vicinity of the Fermi points constitute one Dirac fermion "J" = (tpR, iPl)- 

The marked difference between this, so-called Luttinger, behavior and the Fermi liquid one prompts fundamentally 
important questions pertaining to the possibility of a similar behavior in D > 1 and/or in the presence of long- 
ranged electron-electron interactions. While in the case of the short-ranged interactions the possibility of the D > 1 
Luttinger-like behavior is likely to be limited to the infinitely strong coupling limit, the long-ranged forces appear to 
be capable of destroying the Fermi liquid even at finite couplings. As the best studied example of this kind, the model 
of degenerate non-relativistic massive fermions (T <C /i <C mc 2 ) which are minimally coupled to an abelian gauge field 
was found to have a distinctly non- Fermi liquid behavior Q], although the latter appears to be quite different from 
$H ' the Luttinger one j|. 

More recently, there has been an upsurge of interest in the relativistic counterpart of this model which is a zero- 
density (n = 0) system of the N- flavored relativistic Dirac fermions coupled to an abelian gauge field which is described 
by the standard action of the Quantum Electrodynamics (QED) 

- r N - i 

S[9, % A] = / dx[J2 */(*7/A + 1^ ~ + ^(^A, - d„A,) 2 ] (1) 
/=i 9 

where, for the sake of completeness, we also included a finite fermion mass m. 

Among the previously discussed examples of the 2D condensed matter systems that support the Dirac-like low- 
energy excitations and allow for such an effective description are the so-called flux phase in the planar quantum 
disordered magnets |^,^| and the layered disordered <i-wave superconductors with strong phase fluctuations proposed 
as an explanation of the pseudogap J^|| and insulating (spin density wave) phases of the high-T c cuprates. Also, 
the non-Lorentz- invariant version of QED2+1 was shown to provide a convenient description of the normal semimetalic 
state of highly oriented pyrolitic graphite |§p[ . 

The number of the fermion flavors N depends on the problem in question, although it is not necessarily equal to 
the number of different conical Dirac points in the bare electron dispersion of a lattice system. In all of the previously 
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discussed 2D examples N — 2 is a number of the electron spin components, while the number of the conical 

points turns out to be either two |^|,^| or four which merely forces one to use the four-component Dirac fermions 

and the corresponding (reducible) representation of the 7-matrices 7^ = <r M <g> 03 constructed from the triplet er M of 
the Pauli matrices. 

In the abovementioncd condensed matter-related applications, the effective gauge fields serve as a somewhat exotic, 
yet often more convenient, representation of such bosonic collective excitations as spin or pairing fluctuations, while 
the Dirac fermions correspond to the auxiliary fermionic excitations such as, e.g., spinons "topological" fermions 
, and so forth. Generically, the quantum mechanical amplitudes describing such degrees of freedom turn out to be 
gauge-dependent, while all the physical observables which experimental probes can only couple to must be manifestly 
gauge-invariant . 

Among such gauge-invariant amplitudes, is the one containing a phase factor (sometimes referred to as a "gauge 
connector" or a " parallel transporter" ) 

G? nv (x,y) =< *(x)e l fr A ^( y ) > (2) 

whose suggestive form makes it tempting to identify Eq.(2) with the physical electron Green function (in spite of its 
being gauge- independent, the function Gj nv explicitly depends on the choice of the contour F). 

To this end, it was conjectured B that by analogy with the problem of the compressible Quantum Hall Effect 
described by yet another kind of the 2D auxiliary (this time, non-relativistic) fermionic quasiparticles, the so-called 
composite fermions, interacting with the statistical Chern-Simons field []l2|| , the electron Green function is given by 
Eq.(2) with the contour T chosen as a straight line between the ending points x and y. 

Furthermore, it was argued in Ref. [Q that in the case to = and at energies and momenta which are small as 
compared to the bandwidth and the inverse lattice spacing, respectively, the gauge-invariant amplitude (2) features 
the Luttinger-like behavior with a positive exponent rj (hereafter we use notations qp = q^p^ and p = 7/4V) 

GL(P)-P/P 2 ^ (3) 

which was also invoked in[|J^| to explain the experimental data on angular-resolved photoemission spectra (ARPES) 
in the high-T c cuprates |13 |. 

In the general case of a D-dimensional condensed matter system which possesses a number of isolated Fermi points 
located at fcpi, the conjectured behavior (3) corresponds to the algebraic suppression of the electron propagator at 
long times and distances 

GL(*)~E e ^^^ ( 4 ) 

i 

where the sum is taken over all the Fermi points. 

In the present paper, we employ a functional integral technique to compute the function (3) and discern the true 
nature of its singular behavior (if any) . This approach which had been pioneered by Schwinger and later advanced 
by a number of other authors (see, e.g., [ ]iO||ll| ] and references therein) exploits a functional integral representation 
of the exact solution of the equation for G inv (x, y\A) as a functional of an arbitrary configuration of the gauge field 
A(z). Subsequently, by averaging over the gauge field, one obtains a sum of all the multi-loop diagrams with no 
couplings between the fermion polarization insertions into the gauge field propagators and the open fermion line 
corresponding to the fermion's propagation between the space-time points x and y. Likewise, in the case of a generic 
multi-fcrmion amplitude, the allowed graphs can only contain open fermion lines which connect the incoming and 
outgoing asymptotical fermionic states, provided that the fermion polarization has already been absorbed into the 
gauge field propagator. 

This approach can be viewed as a systematic improvement of the celebrated Bloch-Nordsieck model where all the 
spin-related effects are ignored which makes this model exactly soluble but restricts its applicability to the infrared 
(IR) regime \p 2 — m 2 | -C m 2 near the fermion's mass shell. 

We emphasize that the IR regime can only exist if the fermions are massive, while in the massless case the entire 
region below the upper cutoff A (which is set by the conditions of the applicability of the effective QED-\ike description 
itself) falls into the opposite, ultraviolet (UV), regime which, in the case of a finite fermion mass, is defined as 

\p 2 — TO 2 I ^> TO 2 . 

The rest of the paper is organized as follows. We first describe the Schwinger's functional technique and inves- 
tigate both the IR and UV asymptotics of the ordinary (gauge-dependent) fermion Green function in the general 
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D-dimensional case. Then, after having compared our general formulas with the well known 3D results as well as 
with the partially known 2D ones, we proceed with the gauge-invariant fermion amplitude proposed in Ref. and 
ascertain its true behavior. We conclude our analysis with a discussion of the alternatives to the previously suggested 
form of the physical electron propagator as well as to the fits to the ARPES data [|l3| exploiting the QE D2+i-re\&ted 
scenarios. 



II. FUNCTIONAL INTEGRAL REPRESENTATION OF FERMION AMPLITUDES 



The conventional fermion Green function is given by the (properly normalized) functional integral over the fermion 
and gauge field configurations 

G(x,y) =< V(x)$(y) >= J DfflD[*]D[A]V(x)lt(y) exp(iS[tf,*, A]) (5) 

Upon integrating the fermions out, one arrives at the expression 

G(x,y) = J D[A}G(x,y\A)exp(iS eff [A\) (6) 

where the effective action of the gauge field includes the fermion polarization 

r»i 1 f . , n a r, a , det [id + A — ml 

S effW = ^2 / Md»A v - d v Arf + In — i— , -1 

J det[io — m\ 

= \j d *J dyA^(x)V-i(x - y)A u {y) + ... (7) 

By neglecting all but the gaussian term in (7) one excludes from consideration any processes of "light-light scattering" 
and alike. Thus far, none of the beforementioned effective QED-\ike descriptions of the condensed matter systems 
has gone anywhere beyond this common approximation. 

Nonetheless, the gauge field is not completely quenched, as one still accounts for the quadratic polarization H(q), 
resulting in the gauge field propagator which, in the covariant A-gauge, assumes the form 

= 2 / nr + (A - 1)^) (8) 

q- + U(qj q z 

In turn, the fermion Green function G{x,y\A) computed for a given gauge field configuration obeys the equation 

[id + A(x)-m]G(x,y\A) = S(x-y) (9) 
Its formal solution can be written in the form of a quantum mechanical (i.e., single — particle) path integral ]l0[ | 

G(x,y\A)=-i / dse' s( -- m2+i5) [id + A(x) + m] / D[a]8{x - y - 2 / a(T 2 )rfr 2 ) 
Jo J Jo 

e M-iJ Q dT(a 2 {T)-{2a^(T) + a^id I/ )A fl {x-2 J a(n)dn))] (10) 

where a^ u = [ r y tl ,"f L ,]/2 and 5 — * + . The integral over the fermion's momentum a(s) as a function of the proper time 
s parameterizing its space-time trajectory is normalized in such a way that 

J D[a]exp(-z J a 2 (r))dr = 1 

Next, we perform functional averaging over different gauge field configurations with the use of Eq.(7), then Fourier 
transform Eq.(10) to the momentum representation, and finally switch to the integration over the fluctuating part of 
the total fermion's momentum v(s) = a(s) — p, thus obtaining 

G(p)=-i dse ls{p2 - m2+lS) / L>[v]expH / v 2 (T)dr)[p + m + M(s|v)] exp(i*(a|v)) (11) 
Jo J Jo 
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In this expression, the terms which are odd in A(z) contribute to the gauge invariant (see below) part of the mass 
operator 

M(S|V)= / (o^T^)^^^^ (12) 
while the even ones stem from the exponential of the (gauge-dependent) "phase factor" 

$(s|v) = J ^p5+T ^«'(g) J dri J d,T 2 (2vf_ t (Ti) + 2p ll + (T t _ ta q a ){2v u (T 2 ) + 2p u - Vvpqp) 

2ipq(ri-T 2 )+2i P 1 qv(T 3 )dr 3 . . 

e Jt 2 (13) 

In the above expressions, the integrations over the proper time parameters n are ordered according to the order of 
their appearance in the products of the non-commutative factors {2v^{t{) + 2p^ ± a^q^). 

III. INFRARED BEHAVIOR 

By using Eqs. (11-13) one can readily determine the IR behavior of the fermion Green function. With its momentum 
satisfying the condition \p 2 — m 2 \ <C m 2 a fermion behaves as a heavy particle whose velocity remains essentially 
unchanged atfer emitting and absorbing an arbitrary number of the gauge field quanta. Therefore, the Green function 
receives its main contribution from the fermion trajectories close to the straight-line path (which only coincides with 
the semiclassical trajectory in the case of a time-like separation between the ending points (x — y) 2 > 0). 

This allows one to neglect the fluctuations of the total fermion's momentum with respect to its average value p, in 
which case the mass operator introduces only a small correction 

f da 1 - e 2iqps 1 Id 2 — m 2 l 

M/ * (s|v) = 1 J (^pn^fa)w qp = m— 2 ) ~ P ] - m2 « p (i4) 

In deriving (14) we took into account that a characteristic value of the parameter s ~ \p 2 — m 2 | _1 is determined by 
Eq.(ll) and the fact that the integral (14) receives its main contribution from small transferred momenta q < 1/sp ~ 
\p 2 - m 2 \/p < p. 

In contrast, the integrals over Tj in the gauge-dependent IR phase factor are formally divergent. They must 
be tackled by first computing the momentum integral and then applying the so-called "ribbon" regularization |Tl| l 
p(Ti — t 2 ) — * p(ti — t 2 ) + 1 with (pi) = and |1| = 1/A which yields the expression 

= l9 2 I D [ S dn I" dr 2 [(D-2 + \)— P \ 2(A-1) f (Tl "^i'j 

|p(ri-T 2 )+lr |p(n - t 2 ) + 1| 4 
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= ig 2 I D [{D - 2 + X)(-(s P A) - In(apA)) - 2(A - l)(-(«pA) - ln(spA))] (15) 

In the massive case, the linear divergence of $(s) would be routinely attributed to the renormalization of the bare 
mass m — > m + 0(A). After having separated this linear divergence, we observe that the subleading logarithmic terms 
conspire to give rise to the non-perturbative formula 

G IR (p) = -i(p + m )J~dse^ 2 -^(spA)-^ 2 ~ (p 2_JlTsy-^ ^ 

which, near the mass shell, exhibits the anticipated algebraic behavior (3) with the IR anomalous dimension 

Vl R = 2g 2 I D (\- D) (17) 

where 

I D = [2 D ^ D+1 ^ 2 T({D + l)/2)]- 1 (18) 
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Thus, in the 3D case of the conventional weakly coupled QED^+i we recover the well-known IR exponent (see, e.g., 
@) 



lfg = i^(A-3) (19) 

which vanishes in the so-called Yennie's gauge A = 3 {t]ir is also known to be zero in some non-covariant gauge, such 
as the Coulomb gauge qA — 0). 

In the (parity-even) 2D case which is of a particular interest in view of its condensed matter-related applications 
H the weak coupling regime turns out to be intrinsically unstable against the effects of the fermion polarization. 
In fact, for q < Ng 2 the gauge propagator is totally dominated by the fermion polarization which, for N 3> 1, is given 
by the one-loop term 



(20) 



and the gauge field propagator reads as 



2^( ff ) = _* (^ + (A-1)S^) (21) 

Instead of the bare coupling g, it is l/N that now becomes a parameter of the perturbative expansion. We note that 
above the momentum scale Ng 2 no further logarithmic corrections are generated, so that the latter is now playing 
the role of the UV cutoff. Nonetheless, for the sake of uniformity of our presentation, in the following discussion we 
will continue using the notation A and the label UV for the range of momenta m -C q ^ A = Ng 2 . 

It is also worth mentioning that, owing to the parity conserving structure of the reducible four- fermion representa- 
tion, the radiative corrections generate no Chern-Simons terms. 

Using (21) we obtain a coupling-independent anomalous exponent 

rt° R =^-% (22) 

thus discovering the 2D analogue (A = 2) of the 3D Yennie's gauge. 

Notably, the IR wave function renormalization assumes the anticipated power-law form, in full accord with the 
physical origin of the IR singularity. The latter is known to stem from the processes involving independent emission 
and absorption of an arbitrary large number of soft gauge quanta. Due to their uncorrelated nature, these multiple 
" bremsstrahlung" events obey a Poisson distribution formula, hence the appearance of the factorials in the statistical 
weights, resulting in the natural exponentiation of the lowest order (~ g 2 lnA) correction. 



IV. ULTRAVIOLET BEHAVIOR 



The Schwinger's functional technique is also capable of exploring the UV regime (\p 2 — m 2 \ ^> m 2 ) which is the only 
regime of interest present in the massless case. Despite the fact that the procedure is straightforward, there seems to 
have been no such a systematic attempt made in the past. 

Technically, the UV behavior is more difficult to analyze, because the path integral (11) is no longer saturated 
by the trajectories close to the semiclassical straight line. In fact, the relevant paths can strongly deviate from the 
straight-line one, for they suffer no exponential suppression, unlike in the IR regime. 

Despite the fact that the functional integration over v(s) can no longer be carried out exactly, one can instead 
resort to the formula 

DMe-*f v2dT+F M=e< F > f DMe- l f v2dT jr {F[vh < F>)2 (23) 

n=0 

where < F >= /£>[v]e" i ./ V ' lT F[v]. 

Eq.(23) has been extensively used, e.g., in implementing the Feynman's variational principle in the polaron and 
related problems. Expanding (11) to the first order in T>^ u {q) we obtain 
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SiG uv (p) = -i dse ls{p2+l5) [< M UV (s) > +ip < <S> uv {s) >] (24) 



The functionally averaged mass operator (12) is now determined by the transferred momenta q ^> p ~ and it 

needs to be computed only to the first order in p 

f da 1 — e is (9 2 + 2 qp) D 

< M uv (s) >= i J —^L_ V^(q) q2 + 2qp 7m + 2P» ~ <7»\Q\) = 2g 2 pI D ^-^\n(sA 2 ) + ... (25) 

Notably, Eq.(25) is independent of the gauge parameter. In contrast, the averaged phase factor (13) which can be 
calculated in the p — > limit 



<*w(*) >=/ J^5+T V A<1) 



.1 - e ^(<z 2 +2qp) + is ( q 2 + 2qp ) 



(g M + 2p M + cr AtQ <7a)(g J y + 2p„ - cr^g^) - isS^} = -g I D (D + A) ln(sA ) + ..., (26) 



l (g 2 + 2qp) 2 1 --ff -vfjii-'/ --"imisi 2 

does manifest a dependence on the gauge parameter. Combining (25) and (26) together, we obtain the total correction 
to the Green function 

r°° ■ , i , ( da 1 - e is (q 2 + 2 ip) 

SiG uv ( P ) = J dse"<* +tS} J (^pTT^(g)[ g2 + 2qp b&u + 2Pu °„x<lx) 

I _ e »s(9 2 +2qp) _|_ i s (q 2 -|- 2qp) 

(g 2 + 2qp) 2 ^ + + (7 ^ aqa ^ q " + 2pv ~ cr ^^) ~ ?; P S< W 

, 9 2 p , £>(3-£>) A 2 

= -5-— ^ FiTl Aln — +... (27) 

z yr L> + 1 p A 

By using the identity 

P7m(p + ?)7vP = p(<7a< + 2 ?V + ^iMxQodiqv + 2p v - dvpqp) - JfiP 2 (q v + 2p„ - <r v \q\) - S^pip + q) 2 

and integrating in (27) over the proper time s prior to the momentum integration one can also check that the correction 
given by Eq. (27) exactly reproduces the one- loop result of the conventional diagrammatic expansion 



SiG uv (p) = -1 \ - 4 r ' MP + ®1»P (28) 



dq V^(q) 
(2tt) d +1 p\p + q) 

Instead of expanding Eq.(ll) to higher orders in T>^ u (q) one can perform summation of the leading (g 2 lnA) n terms 
by virtue of the standard renormalization group equation which, reflects the scaling properties of a generic two-point 
amplitude (gauge invariant and non-invariant alike) under the change of the upper cutoff [ jl4| 

[A IK " m Wg + ^)}pGuv(p; A; 9) = (29) 

where the leading order dependence of the anomalous dimension of the fermion Green function on the renormalized 
coupling strength g is given by the explicit form of the first order correction (27) 

d 

v(g) = - a qx p SiGuv (p; a - b )I p =a (30) 

while (3(g) = Adg/dA\ p= \ = 0, and, therefore, the coupling strength retains it bare value g = g, for as long as the 
dynamics of the gauge field is considered quenched. 

The solution of Eq.(29) suggests that the first logarithmic correction (27) merely gets exponentiated, thus yielding 
the algebraic behavior controlled by the UV exponent 

Vuv=9 2 Id(X + ^§^) (31) 
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Further corrections to Eq.(31) require one not only to extract the subleading corrections of order <7 2 ™lnA from the 
n'^-order terms in the expansion of Eq.(ll) in powers of v(s) and account for the improved fermion polarization H(q) 
but also to proceed beyond the quenched approximation (7) for the effective action of the gauge field. 



In the weakly-coupled 3D case Eq.(31) reproduces the well known result 14 

V 3 u D v = (32) 
while in the 2D case it yields the coupling-independent UV exponent 

^ = sAv (3A ~ 2) ' (33) 

in agreement with the result obtained in [fl5| . 

V. GAUGE INVARIANT FERMION AMPLITUDE 

After having tested our formalism against the known examples, we turn to the proposed candidate for the physical 
electron propagator which is given by Eq. (2) with the straight-line contour T 

GL(*>V) = I D[A]G(x,y\A)exp(i f dz^(z)) exp(iS eff [A}) (34) 



v 



Proceeding by analogy with the derivation presented in Section II, one readily obtains Eq.(ll) where Eqs.(12) and 
(13) are replaced, respectively, with 







/o / \ i n \ 2ipqr+2ir/s f" qv(ri)dri f d>T2 , /oe\ 

-(2^(r 2 ) + 2p v - a vX q\)e ™ Jo ^ j j (35) 

o s 



and 



S fTi 



®inv{s\v) = J t^Td+T ^A"^) J dTl j dT 2 (2v f _ t (Ti) + 2p fl + (T tia q a ){2v u (T 2 ) + 2p v ~ o-vpqp) 

, 2ipq(Ti-T 2 )+2i qv(r 3 )dr 3 ^_ 2 T ( ' T4 e 2 'P ( l( T 3- T i)+ 2i ( T 3- T 4)/ s /„* 1v(T 5 )drs 

Jo s Jo s 

f S dr% 2ipq(s-Ti)+2i f" qv(T 4 )eiT4-2ipqT 3 -2iT 3 /s P qv(r 5 )dT 5 
— 2 / e Jt i J " 

Jo s 



(36) 



In the IR regime the path integration can still be carried out exactly by simply neglecting v(s) with respect to the 
average fermion's momentum p. In the same approximation as that used in Section III (which is only justified in the 
vicinity of the mass shell, provided that m^O), one readily obtains 

M mvJR (s) = 2 J —fpLj-lVfo) [ dr ltlPv [e 2 ^ s -^ - e 2 ^] = (37) 

and 

$mv,m(s) = -i J , 2 ^y }+1 1>nv(q)PnPu J drx J dr 2 

j e 2i P q(rx-r 3 ) + 2_ f dT3 f 3 dr 4 e 2lpq(T3 - T4) - - / dr 3 e 2lpclis - Tl - T3) } = 
s Jo Jo s Jo 



(38) 



Thus, as first pointed out by the authors of Refs. p0| , in the IR regime the gauge-invariant propagator (34) retains a 
simple pole 
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p -J- 777 

G m,,/flW « p 2_ m 2 +lS ( 39 ) 



hence, rj inv j R = 0. 

By comparing (39) and (17) one can also deduce the IR anomalous dimension of the exponential factor 
exp(i J dx^A^) itself 

Ve Xp ,iR = '2g 2 I D (D-\) (40) 

which of course vanishes in the Ycnnie's gauge. 

Next, going over to the UV regime and expanding Eqs.(35) and (36) to the first order in T>^ v {q) we arrive at Eq.(24) 
where the functional average of the gauge-dependent phase factor 

/ dc l g 2 



{2tt) d + 1 q 2 + U(q) 



1 - e ^(? 2 +2qp) + i s ( q 2 + 2qp) p 2 q 
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(,= + 2q p,= (-^ ? - M i-f)-i,\ = Ol r > s) <, (41) 

now exhibits neither linear, nor logarithmic divergence as a function of s, unlike in the case of the non-invariant 
amplitude (see Eq.(26)). In turn, the value of the mass operator 

^ at o- / rf q ff 2 l-e^ 2+2 qP) r . ^ . p^qy qp 

< M inVtUV {s) >=2l / — — — [p~q q i h -^cr^gv 

7 (27r)- L '+ 1 + U(q) q 2 + 2qp qp q 2 q 2 

= 2g 2 pI D ^^\n( S A 2 ) + ... (42) 

appears to coincide with Eq.(25). Thus, it is Eq.(42) that solely determines the correction to the gauge-invariant 
Green function 

SiG mv .uv(p) = -i / dse ls( - p2+l5) [< M mv , uv (s) > +ip < ®i„ v ,uv(s) >} = 2g 2 ^I DlT — ln(^) (43) 
Jo p 2 L> + 1 p 2 

The same result can be obtained by working in the axial gauge nA = defined by the vector n = (x — y)/\x — y\. 
In this gauge, the exponential factor in (34) is identically equal to unity, and the first order correction is given by 
Eq.(28) where one has to use the gauge field propagator 

F q 2 + U(qj {nqj 2 (nq) 

Notably, the result (43) obtained with the use of Eq.(28) is independent of the direction of the vector n, for all the 
terms proportional to n(np) cancel out and only those proportional to p remain in the final expression. 

It is worth mentioning that the integrals in Eqs. (41,42) as well as in Eq.(28) with the gauge propagator (44) are 
all plagued with the spurious poles, such as l/(qp) 1,2 . We handle these singular denominators by resorting to the 
exponential integral representation: l/(qn) = — i J dsexp(zs(qn + iS)). Then, after having performed the Lorentz- 
invariant momentum integration, we carry out the remaining integrals over the auxiliary parameter s with the use of 
the "ribbon" regularization |HJ|. This procedure yields the following logarithmic integrals appearing in our calculation 

dq q^ _ H D rfy m( ,A 2 ^ 



{2tt) d + 1 q D - 1 (p + q) 2 (qn) 2 n 2 V 



dq 



q^q v 



(2tt) d + 1 q D - l {p + q) 2 {qnf 



H D 2nn_n v - S^n 2 A 2 

1 ^ In — ), 

2 n 4 p z 



and 



dq 



q^qx 



(2tt) d + 1 q D + r (p + <?) 2 (qn) 2(D + 1) 



( z£ 2 „4 



8 



One can check that the above expressions are fully consistent with the standard "principal value" prescription for the 
spurious poles, whose advanced form is known in the field-theoretical literature as the Leibbrandt-Mandelstam rule 
(see jl6| and references therein). 

Finally, by invoking the renormalization group equation (29) we find that the logarithmic correction (43) tends to 
exponentiate, thereby resulting in the new UV anomalous dimension 

Vinv,uv = -^9 2 Id d ^ l , (45) 

which appears to be negative. 

Surprisingly enough, we were unable to find in the literature any result pertaining to the weakly-coupled 3D abelian 
gauge theory (e.g., conventional QED3+1), in which case Eq.(45) yields 



Vinv,UV — ~g^2 (46) 



Nevertheless, we did find some comfort in comparing (46) with the exponent which had been previously found to 
control the power-law UV behavior of the non-abelian analogue of Eq.(34) in the 5[/(3)-symmetrical case [17j] 

3D,SU(3) (47) 

By construction, Eq.(47) is proportional to the quadratic Casimir operator in the fundamental representation of the 
color group which, in the case of SU (N) , equals 



1 N 1 /V 2 -1 
c F = — V tr(T a T a ) = (48) 

0=1 

Evaluating (48) for 517(3) we obtain c F — 4/3 and, upon separating this factor out, recover the result (46) 
pertinent to the abelian case (with the electric charge e substituted for g). 

Likewise, by using Eq.(21) we obtain the anomalous exponent which controls the gauge- invariant propagator in 
QED2+1 

32 

^nv,UV = (49) 

which is negative, contrary to the result of Ref. J4] and in agreement with the sign (albeit not the magnitude) of 
the exponent quoted without derivation in Ref. fel. However, it remains to be seen whether the exponentiation of 
M^{s) as well as vanishing of 3>^(s) still hold beyond the leading 1/N order. 

Lastly, by comparing Eqs.(31) and (45) one can also deduce the UV anomalous dimension of the exponential factor 
exp(i / A^dz^) 

Vex P ,uv = -9 2 Id(D + A) (50) 

Interestingly enough, for A = —D this exponent equals zero, and the UV anomalous dimension of the non-invariant 
propagator coincides with (45), in agreement with the observation made in the 3D non-abelian case ffffl. 



VI. DISCUSSION 



Our calculation demonstrates that in the massless case the gauge invariant Green function (34) appears to decay 
slower than the bare one, in a marked contrast with the previously conjectured Luttinger-like behavior. In this 
concluding Section, we make an attempt to rationalize these findings, although we refrain from making any final 
judgement on their physical implications. 

Albeit somewhat counterintuitive, the found UV behavior is not totally incomprehensible. In fact, the generic 
behavior of an invariant fermion amplitude is manifested by the asymptotic formula 

G\ nv {x) ^eM-C\x\K + n\n{\x\K)) (51) 
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where C > 0, and the expression (51) decays with \x\ exponentially, regardless of the sign of 77, because the logarithmic 
term in the exponent is subleading to the linear one. However, in a renormalizable gauge theory where the gauge 
invariance is reinforced throughout the whole process of renormalization, the latter would be routinely cancelled out 
by counterterms, which leaves behind the logarithmic term of (potentially) either sign. 

This situation would change, however, should one choose to relaxe the condition of renormalizability at the expense 
of the gauge invariance, since the radiative corrections to the action (1) generically produce a finite mass of the vector 
field A^. Loosely speaking, the situation would then resemble that in the Schwinger's QED\ + i where the gauge field 
acquires a mass M ~ g, and the analogue of Eq.(34) behaves as 

G[ nv (x) ~ exp(-i[ln(Mx) + K (Mx)] - Mx) (52) 

It is worth noting that, should one decide to intentionally disregard the exponential factor e~ Mx , Eq.(52) would 
appear to exhibit a power-law decay ~ 1/s/x at x ^ 1/M, thus siggesting -q}^ = —1/2. 

We mention, in passing, that the exponential, rather than a power-law, behavior has also been found in the 
problem of Dirac fermions in the presence of a static random vector potential (A(x) = (0, a(f))) which allows for an 
asymptotically exact solution in the ballistic regime of large fermion energies [p_8| . 

Conceivably, in some of the abovementioned physical applications of QED 2 +i with N — 2 the problem of the slow 
space-time decay of the gauge invariant amplitude (34) can be thwarted by a spontaneous development of a finite 
fermion mass, in which case the behavior of G inv (x) at large x will be governed by the (free) IR asymptotic (39) 
instead of the UV one. However, the intrinsic propensity of the 2D Dirac fermions in QED2+1 towards generating a 
finite mass (usually referred to as the phenomenon of chiral symmetry breaking) is believed to occur only at sufficiently 
small N < N c Jl9| . While in the case of the Lorentz-invariant action (1) the critical number of flavors N c was found 
to be as low as 3/2 |L9|], the Lorentz-(or even rotationally-) non-invariant generalizations of the action (1) are still 
awaiting to be fully explored. 

To this end, the authors of Refs. Q conjectured that the critical value N c in the QED-like description of the quantum 
disordered planar d-wave superconductor may become greater than two due to the lack of rotational invariance. On 
the other hand, in the finite-temperature counterpart of the 2D chiral symmetry breaking transition in the (spatially) 
rotationally-invariant effective theory of a single layer of graphite N c was found to be further reduced as compared to 
the Lorentz-invariant case Q. 

However, should one insist on maintaining both the gauge and Lorentz invariances of the renormalized gauge field 
action, the problem of the slow spatial decay of the alleged physical electron propagator (34) associated with its 
negative UV anomalous dimension (45) could not be resolved without re-examining the "minimal" form of this Green 
function. In fact, the task of constructing the proper gauge transformation which converts the auxiliary Dirac fermions 
into the physical electrons may not be limited to a particular choice of the contour T in Eq.(2) but may also require 
one to modify the phase factor itself. 

It is worth noting that in the previous calculations of the "zero-bias anomaly" in the tunneling density of states 
in the compressible Quantum Hall effect [Q , the construction of the electron Green function, albeit seemingly given 
by the same Eq.(34) with the contour T now chosen along the temporal axis, was, in fact, more involved. Indeed, in 
the semiclassical approximation employed in [p"2|| , the gauge field dependence of the exponential factor exp(i J dz^A^) 
would have been exactly compensated by that of the non-gauge invariant Green function G(t, 0|A), thus making the 
functional average of the product of the two behave essentially as in the absence of any gauge coupling. 

Nevertheless, the electron density of states computed in jL2j appeared to be strongly affected by the Chern-Simons 
gauge fluctuations which can be traced back to the fact that, in addition to the abovementioned factors, the electron 
Green function happened to contain yet another factor, the exponent of the saddle-point value of the effective action 
of the Chern-Simons gauge field. It was, in fact, this factor which was solely responsible for the strong suppression of 
the tunneling density of states, consistent with the physical interpretation of the Chern-Simons field as representing 
the effect of the Coulomb coupling in the presence of strong magnetic field. In light of the fact that in the problem 
at hand the time reversal symmetry remains unbroken, no such an additional factor can be readily incorporated into 
the naive form of the electron propagator (34) . 

In order to further elaborate on this point, we mention yet another example demonstrating the sensitivity of a 
generic gauge invariant amplitude to the details of its construction. To this end, we recall the original Dirac's idea of 
explicitly constructing a "dressed charge" corresponding to a physical electron by means of the gauge transformation 

$ P hys(x) = exp(i / dyx^ix - y)A fl (y))V(x) (53) 
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where the vector function Xn( x ) obeys the equation d^x^ix) = 5{x). In the time-independent Shrodinger operator rep- 
resentation, the originally proposed transformation from the bare fermions to the physical electrons was implemented 
as a space-like Dirac string between the location of the fermion and an infinitely remote point 

Xo = 0, Xi =< x\^d l \y> 

The Fourier transform of the electron propagator < ^ p hys(%)^phys(y) > is IR finite (see Eq.(39)) and undergoes 
multiplicative UV renormalization at a single point p^ — (m, 0) on the mass shell corresponding to a static charge 
p6[ , in agreement with the general expectation that the absence of any singularity other than a simple pole is 
characteristic of the propagator of an exact eigenstate with the quantum numbers of electron. 

It was shown in that in the case of a dressed charge moving with a finite velocity u the above phase factor 
needs to be further modified 

^ P hys(x\u) = exp(i 7 J d D - l y^dy\\ < x\^\y > [7 _2 <9||^|| + d ± Aj_ - uE])*(x) (54) 



where 7 = — ip and both the parallel and perpendicular components of A are determined with respect to 

the velocity vector. As shown in [fl6|| , Eq.(54) gives rise to the operator whose propagator is both IR-finite and 
UV-renormalizable at p^ = 7717(1, u). 

Such a strong dependence on the exact details of the construction of the phase factor appearing in the gauge 
transformation (53) indicates that the true electron Green function may well be quite different from Eq.(34). In 
particular, it remains to be seen whether one can at all find an alternate form G\ nv {x) which would decay faster than 
the bare propagator. Given the intellectual appeal of the QED2+1 picture, such an endevour is definitely worth the 
effort, and a further investigation into this possibility is currently under way. 

Should, however, the sought-after Luttinger-like behavior fail to occur even in the modified prototype of the electron 
propagator, one can still consider an alternative approach to the quantum disordered d-wave systems, e.g., the one that 
was put forward in the context of the scenario of a second pairing transition in the 2D superconducting phase |2(]] . In 
pl| , apart from fully idenifying the true nature of this transition and its critical properties (the specific predictions of 
Rcf. |2l[] for the critical exponents are roughly consistent with the recent tunneling data in Ca-doped YBaCuO p2|), 
it was further speculated that it might be possible to extend the effective Higgs- Yukawa theory of the nodal fermion 
excitations coupled to the fluctuations of the secondary order parameter of either id xy or is symmetry well into the 
pseudogap phase. Rather than a global superconducting coherence, this would only require the presence of a local 
parent g?^ ^2 -wave order. If this speculation prove valid, it can provide a viable alternative to the QE D2+i-ba,sed fits 
to the ARPES data , since in the Higgs- Yukawa theory the anomalous dimension of the Dirac fermions is indeed 
positive HHH. 

To summarize, in the present paper we applied the Schwingcr's functional integral representation of the fermion 
amplitudes to the analysis of both the infrared and ultraviolet asymptotics of the conventional (non-gauge invariant) 
fermion Green function and a particular gauge-invariant amplitude (34) . 

In the IR regime, this method provides a substantial improvement with regard to the spinlcss Bloch-Nordsieck 
model or the customary semiclassical (eikonal) approximation, since it preserves the exact spinor structure of the 
fermion amplitudes. Moreover, the intrinsic "exponential" form of the Schwinger's integral representation facilitates 
truly non-perturbative calculations. 

In the opposite, UV, regime, the method allows one to naturally separate between the gauge invariant and non- 
invariant contributions to the mass operator and systematically compute the higher order contributions into both kinds 
of terms. For a specific class of problems, including the amplitudes given by Eq.(2), it has a significant advantage as 
compared to the conventional diagrammatic technique which is not particularly well suited for such calculations, for 
the very rules of the diagrammatic expansion turn out to be amplitude-specific and depend on a particular choice of 
the contour V @. 

To our surprise, the previously suggested "minimal" form of the physical electron Green function (34) was found 
to manifest a negative anomalous dimension, contrary to the much-anticipated Luttinger liquid behavior. The impli- 
cations of this observation were discussed, some of them pertaining to the applicability of Eq.(34) and other to the 
possible alternatives to the QED2+i-like description of the ARPES data in the high-T c cuprates. 
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